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Abstract. In a recent paper [A. Santos, G. M. Kremer, and V. Garzo, Prog. Theor. Phys. Suppl. 184, 31-48 (2010)] the 
collisional energy production rates associated with the translational and rotational granular temperatures in a granular fluid 
mixture of inelastic rough hard spheres have been derived. In the present paper the energy production rates are explicitly 
decomposed into equipartition rates (tending to make all the temperatures equal) plus genuine cooling rates (reflecting 
>0 , the collisional dissipation of energy). Next the homogeneous free cooling state of a binary mixture is analyzed, with 

Cn ■ special emphasis on the quasi-smooth limit. A previously reported singular behavior (according to which a vanishingly 

small amount of roughness has a finite effect, with respect to the perfectly smooth case, on the asymptotic long-time 

^ translational/translational temperature ratio) is further elaborated. Moreover, the study of the time evolution of the temperature 

^4— ( , ratios shows that this dramatic influence of roughness already appears in the transient regime for times comparable to the 

O ■ relaxation time of perfectly smooth spheres. 
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Tj INTRODUCTION 



The simplest model of a granular fluid consists of a system of identical, inelastic smooth hard spheres with a 
constant coefficient of normal restitution [1]. Obviously, the model can be made closer to reality by introducing 
extra ingredients. In particular, polydispersity and roughness are especially relevant because they unveil an inherent 
breakdown of energy equipartition in granular fluids, even in homogeneous and isotropic states [2, 3, 4, 5, 6, 7]. 
In this work a fluid mixture of inelastic rough hard spheres characterized by mutual coefficients of normal restitution 
f— ^ ' {(^ij} ^nd tangential restitution {jitj} is considered. First, a recent derivation by kinetic-theory arguments [8] of the 

[~^ . collisional energy production rates ((^/^ and ^/j") associated with the translational and rotational temperatures (7]"^ and 

^^ I 7:™') is recalled and the energy production rates are decomposed into equipartition and cooling rates. Next, the results 

[^^ . are applied to a binary mixture and the asymptotic long-time values and the time evolution of the three independent 

f^ ' temperature ratios T2 /Tf, T™^ /T™^, and T^'/Jj"^ and are analyzed in the homogeneous and isotropic free cooling 

^D . state Special attention is paid to a paradoxical effect: the asymptotic long-time values of the temperature ratio T2 /T^ 

in the nearly-smooth limit (/3y — >^ — 1) differ from those obtained for purely smooth spheres ()3,, = — 1). 



COLLISIONAL ENERGY PRODUCTION RATES 



Let us consider a multi-component granular gas made of hard spheres of masses {m,}, diameters {cJ,}, and moments 
of inertia {/,■ = ^miafKi}. The value of the dimensionless parameter K", depends on the mass distribution within the 
sphere and runs from the extreme values fc, = (mass concentrated on the center) to k", = | (mass concentrated on the 
surface). Collisions between a sphere of component / and a sphere of component j are characterized by a coefficient 
of normal restitution a,y and a coefficient of tangential restitution j3,y. The former coefficient ranges from a,y = 
(perfectly inelastic particles) to a,y = 1 (perfectly elastic particles), while the latter runs from j3,; = — 1 (perfectly 
smooth particles) to j3,y = 1 (perfectly rough particles). The total (translational plus rotational) kinetic energy is 
dissipated upon collisions unless a,-, = 1 and jS/y = ±1. The particular case of one-component systems has been 
widely studied in the literature (see, for instance, Refs. [5, 6, 7]). 

Let /,(r,-,v,-,a),;f) and fij{rj,\j,a>f,rj,\j,G)j;t) be the one-body and two-body distribution functions, respectively, 
where v, is the velocity of the center of mass and fi), is the angular velocity. By starting from the Liouville equation 
and foUowing standard steps, one can derive the Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) hierarchy. The 



first equation of the hierarchy reads [8] 

dtM\i,mi)+\i-Vfi{yi,mi) = Y.M'fi:0>i\fij], (1) 
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Here, 7,,- is the colHsion operator, a,, = (a, + <7,)/2, v,, = v, — v/, rj*^ = r, ± CJ,j5, and the double primes denote pre- 
colHsional quantities giving rise to unprimed quantities as post-colhsional values. The explicit form of the restituting 
collision rules can be found in Ref. [8]. 

Given an arbitrary one-body function V^,(v,, O,), we define its average as 

(V/,(vi, <»,•)> = - / d^i j dOii WiiM, oai)fi{\i, (»,•), rii = / d\i / dO)ifii\i, (B,), (3) 

where n, is the number density of component /. The rate of change of the quantity V^,(v,, O,) due to collisions with 
particles of component j is 

^ij[\l/i{vi,a>i)] = - fdvi I da>i\i/i{\i,a>i)Jij[\i,oii\fij]. (4) 

In particular, the partial temperatures associated with the translational and rotational degrees of freedom are 

Tr^'^{{y,-uf), rr^^ico^), (5) 

where u = Y.if^if^ii'^i) /Hi'^if^i is the flow velocity. The corresponding collisional rates of change define the (partial) 
coUisional energy production rates as 

j j i i 

The global temperature of the gas and its associated net cooling rate are 

where n = ^, n, is the total number density. Note that, as usually done in the literature on granular gases, the Boltzmann 
constant has been absorbed in the definitions of Tj"^, Tj™', and T, so that these quantities have dimensions of energy. 

The energy production rates (^f- and <^,''?' defined in Eq. (6) do not have a definite sign. They can be decomposed 
into two classes of terms: equipartition rates and cooling rates (see Fig. 1). The equipartition terms, which exist even 
when energy is conserved by collisions (a,; = 1 and j3,/ = ±1), tend to make temperatures equal [10]. Therefore, they 
can be positive or negative depending on the signs of the differences Tj'' — T? Tj™' — T™', and 7j'' — Tj™'. On the other 
hand, the cooling terms reflect the collisional energy dissipation and thus they are positive if a,-, < 1 and/or |j3,,| < 1, 
vanishing otherwise. Only the cooling terms in (^?J and <^'?' contribute to the net cooling rate ^. 

The quantities ^,j[v/,(v,, ffl,)] are functionals of the two-body distribution function fij and therefore they are 
generally rather intricate. Now, let us imagine that, instead of the full knowledge of /],-, we only have at our disposal 
the local values of the two densities (n, and nj) and the four partial temperatures (7j"^, 7^™', T? and T™'). Then, we 
can get reasonable estimates of ^,,[v/,(v,, (0,)] from the replacement 
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where Xij is the contact value of the pair correlation function. Equation (8) can be justified by maximum-entropy 
arguments, except that here a Maxwellian form for the rotational probability densities /^™'((B,) and /y°'((By) is not 
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FIGURE 1. Scheme on the two classes of contributions (equipartition rates and cooling rates) to the energy production rates (^P- 
and ^Jf^ characterizing the effect on 7^"^ and 7^™', respectively, of collisions with particles of component j [9]. 



needed. The expressions for ^/J and (^/^' derived from the replacement (8) and the assumption (v,) = (v^) can be found 
in Ref. [8]. Here, those expressions are recast into sums of equipartition and cooling rates as follows: 
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In these equations, Kij = KiKj{mi + nij) / {Kinii + Kjnij), 
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is an effective collision frequency, and we have assumed (o,) = {a)j) ~ 0. The quantities ^^ '"'' represent equipar- 
tition rates. They do not have a definite sign and vanish if all the temperatures are equal. The equipartition rate ^i- 
is always present (even for perfectly elastic spheres, a,y = 1) and tends to equilibrate both translational temperatures. 

The rates ^^ and ^>j do not contribute in the case of smooth spheres (jS/y — —1). The former tends to equilibrate 
the translational and rotational temperatures of component /, while the latter tends to equilibrate both rotational tem- 
peratures but is also affected by the other temperature differences. The quantities i^,'' and Q'^^, on the other hand, are 
positive definite and represent cooling rates. The former only vanishes if the spheres are elastic, whilst the latter only 
vanishes if the spheres are either perfectly smooth (j3,/ = —1) or perfectly rough (j3,y = 1). It is straightforward to check 
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expected, the equipartition rates do not contribute to the net cooling rate. The latter is given by 
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FIGURE 2. Schematic representation of the energy transfer channels in the free cooling state of a binary mixture. The left 
diagram corresponds to perfectly smooth particles (j8 = — 1), in which case both rotational temperatures are completely decoupled 
each other and also from the translational temperatures, so they do not change with time. The right diagram corresponds to quasi- 
smooth particles (/3 > — 1), in which case weak channels of energy transfer exist and also the rotational temperatures are subject to 
a weak cooling. 

HOMOGENEOUS FREE COOLING STATE 

In the homogeneous and isotropic free coohng state the evolution equations for the total and partial temperatures are 



T^tr T^rot 



(16) 



The first equation describes the monotonic decrease of temperature (unless a,y ~ 1 and j3,-; = ±1 for all the pairs). The 
remaining equations show that the temperature ratios Tj"/?^ ^"d 7j™'/r decrease or increase, depending on the sign 
of the differences ^/"^ — ^ and ^ "' — ^, respectively. After a certain transient stage, those temperature ratios eventually 
reach constant values, while the total temperature keeps decreasing following Haff's law T{t) = r(fo) [1 + ^{to)t/2]^ . 
In this asymptotic regime one has equal production rates, i.e., ^f^ — £,2 — ■ ■ ■ — 4"^ — i^[°' — £,2°^ — ■ ■ ■ — ^™', where s is 
the number of components. Use of expressions (9)-( 13) allows one to get an algebraic set of equations whose solution 
gives the asymptotic long-time values of the temperature ratios. 

In the particular case of a binary mixture (s — 2), it is convenient to choose the three relevant temperature ratios as 
T2 jT^, r2™'/r™', and T™^ /Tf . The parameter space is twelve-dimensional: the three coefficients of normal restitution 
(an, a\2, 01,22), the three coefficients of tangential restitution (j3n, j3i2, fe), the two moment-of-inertia parameters (K\ 
and K2), the mass ratio (m2/mi), the size ratio (02/ '^\), the concentration ratio (112/111) and the total packing fraction 
(p =^ (7r/6)(ni(7j'+«2cr|). 

For the sake of illustration, here we particularize to an equimolar mixture where all the spheres are uniformly solid 
and are made of the same material, the size of the spheres of one component being twice that of the other component. 
More specifically, ni ~ n2, an — an — a22 = a, j3n = j3i2 = fe = P, ^i — ^2 ^ 5. O'2/o'i = 2, and m2/m\ — 8. 
Moreover, we consider a dilute granular gas (0 <C 1), so that Xij ~ 1- Thus only the parameters a and j3 remain free. 

In Ref. [8] it was observed that the asymptotic value of the translational/translational temperature ratio T^'' /T^ 
exhibits a peculiar behavior in the smooth-sphere limit j3 — > — 1: it tends to a finite value different from (in fact 
higher than) the value directly obtained in the case of perfectly smooth spheres [2]. Thus, a tiny amount of roughness 
has dramatic effects on the temperature ratio Tj'/T", producing an enhancement of non-equipartition. This singular 
behavior of the case /3 = — 1 in the free cooling state is further elaborated in this paper 

The physical origin of the singularity is the following one. If one strictly has j3 = — 1, the rotational degrees of free- 
dom are "frozen", so that the rotational temperatures are finite while the translational temperatures keep decreasing. 
The rotational/translational temperature ratios diverge but, since the spheres are perfectly smooth, there is no mecha- 
nism transferring energy from the rotational to the translational degrees of freedom. This situation is represented by the 
left diagram of Fig. 2. On the other hand, if /3 = — 1 + e, where < £ ^ 1 , the huge rotational/translational temperature 
ratios are eventually able to "feed" the weak energy channels connecting the rotational and translational temperatures 
(see the right diagram of Fig. 2), thus producing a non-negligible effect ("ghost" effect) on the ratio T^'' /T^ . 

Let us see this in more detail. In the limit j3 — > — 1 one can see that the rotational/translational temperature 
ratios diverge as T^'/r" « d\e^^ and Tj^'/rj" « 02£^- It can be found that the energy production rates behave as 
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where F, X are explicit functions of the indicated arguments. Since both rotational production rates vanish in the 
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FIGURE 3. Plot of the asymptotic long-time values of the temperature ratios Tj /T^ (left panel), ( 1 + j8 ) r|°yrj °' (middle panel), 
and (1 +j8)ri'"°7r/' (right panel) as functions of a for /3 = -0.8, -0.85, -0.9, -0.95, and in the limit /3 -^ -1. In the left panel 
the translational/translational ratio for the strict case of perfectly smooth spheres (j3 = — 1) is also plotted. 
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FIGURE 4. Time evolution of the temperature ratios T^ /Tl' (left panel), T^°yTl°^ (middle panel), and Tl°^ /Tl' (right panel) for 
a = 0.5 and /3 = -0.95, -0.98, -0.99, and -1. The initial state is an equipartition one, i.e., r/''(0) = r2""(0) = r™'(0) = r2™'(0). 



limit e -^- 0, the condition of equal production rates implies that F{\t:^ / 71^02) = 0, f:^{t:^ /T^ .Q^) = 0, and 
Fl°\Tl'/T{', 02/01 ) = FI°\t:^' /Tl'). The first two conditions yield ^^ /T{' and 02- Once obtained, the third condition 
gives 01. In the pure smooth-sphere problem, only the translational degrees of freedom matter and therefore the 
temperature ratio T2 /T^^ is obtained from the condition t,^ ~ £,2 with 1 +/3 = 0, which becomes Fy{T^''/T^,Q) = 
F^'iT^'/Tl'^O). Obviously, the obtained solution differs from that given by F^'iT^'/Tl', 62) = F^\T^'/Tl', §2) = 0. 

In the preceding discussion we have implicitly assumed that the spheres are inelastic (a < 1). In the perfectly elastic 
case (a ~ 1), the scalings for the temperature ratios are quite different and turn out to be T^^ /T^ = 1 + ye, Tl°^ /T^ = 
i?ie, and rj^Vr/'' = ^e. The production rates behave as ^{'/vn ~ '^'Ut)'^^ ^i/vn ~ 'I'"' ("7)e. ^r/vii « <P\°'{^i)e, 
and ^2" 



'/ Vii « <I>2"(i32)e, where 'Pi'^ are linear functions. The coefficient 7 is simply obtained from 'I>['(7) — 't'j (7). 



Once 7 is known, &i and §2 follow from 4>j°'(iJi) = 4>j''(7) and ^"'('^i) = *^9 (7)' respectively. 

Figure 3 shows the quantities T2 /T^, (1 +j3)r2°Yr™', and (1 + fi)Tl°^ /Tf as functions of a for several values 
of /3, including the limit j3 — > — 1. It is clearly observed that the translational/translational temperature ratio in the 
limit j3 ^- — 1 differs from the values obtained in the perfectly smooth case [2]. In the latter case, moreover, the 
rotational temperatures T^^ and T2°^ maintain their initial values and thus the scaled ratio (1 + j3)72°'/r™' vanishes, 
while {\+^)Tl"^/Tl' is indetermined. 

Since the scaled energy production rates lim,^oo(^i 2° /vn go to zero in the quasi-smooth limit /3 -^ —1, it 
is reasonable to wonder whether the paradoxical singular effect, lim,^„ T2 /Tf\n^_-^ ^ lim,^«, ^2'V^i'%=-i' is 
apparent only for times much longer than the characteristic relaxation time of T2 /T^\o^_y In order to investigate 
this possible scenario, the set of three coupled equations d^ \n{Tl' /T{') = - {B,^ - -^f)/ Vi 1 , d^ \n{Tl°^lT{') = - (<^[°' - 
^f)/vii' and dT:\n{T2°^ /Tl^) — — (^9°' — <^i*)/vii has been numerically solved for several values of j3 > —1, the 



solution being compared with that of dr ln{T2/Tl^) = - {^2 ^ ^i)/^i 1 at j3 = — 1 . Here, T = /q dt' Vi 1 (f') is a measure 
of time in units of the accumulated number ofl-lcollisions. Figure 4 shows the time evolution of the three temperature 
ratios for a — 0.5 and some representative values of j3. We can observe that in the perfectly smooth case (/3 = — 1) 
the translational/translational temperature ratio monotonically increases until reaching a plateau value at T « 4, the 
rotational/rotational ratio keeps its initial value, and the rotational/translational ratio increases exponentially with T 
(Haff 's law). The important feature in the cases /3 7^ — 1 is that, while the relaxation time needed to reach the asymptotic 
long-time values actually increases as 1 + j3 decreases, the curves for T2 /Ty and Tl°^ /Ty markedly differ from the 
case of perfectly smooth spheres just after T ss 4. Therefore, the dramatic difference between the quasi-smooth limit 
(/3 —7- —1) and the perfectly smooth case (/3 = —1) manifests itself not only in the asymptotic long-time regime but 
also in the transient regime for times comparable to the relaxation time of perfectly smooth spheres. 

CONCLUSION 

In this paper the homogeneous free cooling state of a binary granular gas made of inelastic rough hard spheres has 
been analyzed. Special attention has been paid to the quasi-smooth limit (j3 -^ —1), where an interesting singular 
phenomenon appears. It turns out that a vanishingly small amount of roughness has a significant effect on the 
translational/translational temperature ratio, with respect to the strict perfectly smooth case (j3 = — 1). This paradoxical 
("ghost") phenomenon is essentially due to the fact that, even if j3 is close to —1, the transfer of energy from the 
rotational to the translational degrees of freedom eventually becomes activated when the rotational temperatures are 
much larger than the translational. 

It is important to remark that this singular behavior of the smooth case is closely tied to the non-stationary character 
of the free cooling state and thus disappears in homogeneous steady states. For instance, in a granular gas driven 
by a white-noise thermostat, the stationarity conditions are {Ty /mi)E,y — {T2 /m2)^2 ^^'^ ^\°^ ~ *32°' ^ ^- ^^ ^^ 
straightforward to check that in this case the ratio T2 /Ty in the limit j8 ^ — 1 coincides with that of j3 = — 1 [4]. 
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